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Recent experiments on Coulomb drag in the quantum Hall regime have yielded a number of sur- 
(~| I prises. The most striking observations are that the Coulomb drag can become negative in high 

^ , Landau levels and that its temperature dependence is non-monotonous. We develop a systematic 

diagrammatic theory of Coulomb drag in strong magnetic fields explaining these puzzling exper- 
QQ ' iments. The theory is applicable both in the diffusive and the ballistic regimes; we focus on the 

experimentally relevant ballistic regime (interlayer distance a smaller than the cyclotron radius Re)- 
It is shown that the drag at strong magnetic fields is an interplay of two contributions arising from 
different sources of particle-hole asymmetry, namely the curvature of the zero-field electron disper- 
sion and the particle-hole asymmetry associated with Landau quantization. The former contribution 
is positive and governs the high-temperature increase in the drag resistivity. On the other hand, 
0^ ' the latter one, which is dominant at low T, has an oscillatory sign (depending on the difference in 

D I filling factors of the two layers) and gives rise to a sharp peak in the temperature dependence at T 

of the order of the Landau level width. 
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PACS numbers: 73.63.-b, 72.10.-d, 73.23.-b, 73.43.-f 



INTRODUCTION 



O ' Coulomb drag between parallel two-dimensional electron systemg^'^ has developed into a powerful probe of quantum- 

, ^, I Hall systems f^i^i^iSiLSiSiifl providing information which is complementary to conventional transport measurements. The 

drag signal is the voltage V developing in the open-circuit passive layer when a current / is applied in the active 

T—i , layer. The drag resistance (also known as transresistance) is then defined by Rd = V/I. As a function of interlayer 

^ ' spacing a, the interlayer coupling changes from weak at large spacings where it can be treated in perturbation theory, 

^^ to strong at small spacings where it can result in states with strong interlayer correlationa^iS. In the present paper 

we will be concerned with the regime of weak interlayer interaction. 

In a simple picture of Coulomb drag, the carriers of the active layer transfer momentum to the carriers of the 

^^ . passive layer by interlayer electron-electron scattering. Due to the open-circuit setup, a voltage V develops in the 

^^ ■ passive layer, which balances this momentum transfer. The phase space for interlayer scattering is proportional to 

C^ . the temperature T in either layer predicting a monotonous temperature dependence Rd oc T^ of the drag resistance. 

^i ■ Moreover, the signs of the voltages in active and passive layer are expected to be opposite (the same) for carriers of 

2 ',. equal (opposite) charge in the two layers>ii It is conventional to refer to the sign resulting for like (unlike) charges as 

^ ' positive (negative) drag. It is worth emphasizing that, as the above considerations imply, the non-zero value of drag 

l^ in the regime of weak interlayer interaction is entirely due to the violation of the particle-hole symmetry. 

Remarkably, experiments show that Coulomb drag behaves very differently from these simple expectations when 
Q ' a perpendicular magnetic field B is applied such that the Fermi energy Ep is in a high Landau level, Ep /fwoc 3> 1. 
^ , {bJc is the cyclotron frequency.) Several experiments^'^' in the regime of weak interlayer coupling observed negative 
ILJ ■ drag when the filling factors in the two layers are different. A more recent experiment! 9- also reveals a non-monotonic 
. ^H , dependence on temperature. While the drag resistivity shows a quadratic temperature dependence at sufficiently high 
j^ ■ temperatures, where drag is always positive, an additional peak develops at low temperatures which can have both a 
J-j positive or a negative sign depending on the filling-factor difference between the two layers. 



Early theoretical work'-^ on Coulomb drag in a magnetic field in the limit of high Landau levels showed that the 
magnetic field may strongly enhance the Coulomb drag, as indeed observed experimentally. On the other hand, the 
calculation of Ref. \VA as well as of a later paper ji^ results in a strictly positive transresistivity, in contradiction 
with the oscillatory sign found in recent experiments. As we discuss in detail below, a general formula for the drag 
resistivity obtained in Ref. [1^ which looks like a natural generalization of the zeio-B resuhiii*i^*i& and also served 
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as a starting point for Ref. [ij, misses an important contribution. This strongly restricts the range of vahdity of the 
results of Refs. Il2ll3 . making them inapplicable under typical experimental conditions. More recent worki^ showed 
that Landau-level quantization can lead to sign changes in drag. However, the results obtained in Ref. [Ij suggested 
that unlike the experimental observation, negative drag should be observed for equal filling factors in the two layers. 
The temperature dependence of the drag resistivity was not studied in 14. 

In this paper, we present a systematic study of Coulomb drag in the limit of high Landau levels. We focus on the 
experimentally relevant limit of well-separated Landau levels (LLs) in which the LL broadening A is small compared 
to the LL spacing htOc- Our starting point is the diagrammatic Kubo formulation of Coulomb dragi^^ for weak 
interlayer interaction. Disorder is included at the level of the self-consistent Born approximationii (SCBA) which 
becomes exact in the limit of high Landau levelsi^. 

Our results are in good agreement with the experimental observations. We find that at high temperatures, the 
leading contribution to Coulomb drag is due to the breaking of particle-hole symmetry by the quadratic dispersion of 
the electrons. This contribution which is analogous to the conventional contribution to drag discussed above, always 
has a positive sign and depends on temperature as T^. At temperatures fcsT <C A, wc find that the dominant 
contribution arises from the breaking of particle-hole symmetry due to the Landau-levcl structure. This contribution 
gives rise to a peak in the temperature dependence and can take on both positive and negative signs, depending on 
the filling-factor difference of the two layers. In particular, the sign is negative for equal filling factors in the diffusive 
regime where the interlayer distance a is larger than the cyclotron radius Re, as was found in Ref. Il4l We find, 
however, that this sign becomes negative in the experimentally relevant ballistic regime (a small compared to Re), in 
agreement with experiment. 

This paper is organized as follows. Sec. ^briefly summarizes the pertinent background on the Kubo approach 
to Coulomb drag as well as on the self-consistent Born approximation. In Sec. IIIII we present the diagrammatic 
calculation of the triangle vertex entering the expression for the drag conductivity, for well-separated LLs, both in the 
diffusive and in the ballistic regime of momenta. In Sec. IIVI we collect the relevant results for the screened interlayer 
interaction. These building blocks are used in SecElto compute the drag resistivity. In this section, we also compare 
our results with experiment. Finally, Sec. IVII contains a summary of our results and a discussion of prospects for 
future research. In what follows, we set h = ks = I- 

II. BACKGROUND 

A. Drag 

Our considerations are based on the Kubo approach to Coulomb dragi^ii^ which expresses the drag conductivity 
crj^(Q, fl) in terms of a current-current correlation function. 



1 r°° I 

^g(Q,^) = ^y^ dte'-^' ([jr^^(Q,i),jf (Q,o)] 



(1) 



where i,j label the components of the drag conductivity tensor, Q, J7 denote the wave vector and frequency of the 
applied field, S is the area of the sample, and j\ denotes the ith component of the current operator in the ^th layer. 
The dc drag conductivity follows by taking the limit 

a,^ = a,^(Q = 0,f]-.0). (2) 

When computing the retarded correlation function appearing in Eq. Q| within the Matsubara technique, the 
leading diagrams in the limit of weak (screened) interlayer interaction C/(q, w) are shown in Fig.^ Analytically, these 
diagrams are given by the expression 

cr^jiiflk) = y^ Tf\c[,iuJn + i^k,ii^n)'rj\q,iuJn,iuJn + i^k)U{c[,iuJn + i^k)U{c[,iLUn)- (3) 

Here, u;„ and ilk denote bosonic Matsubara frequencies and the vector F*-" (q, iujn, ioJm) is the triangle vertex of layer 
I as defined by the diagrams in Fig.|21 Neglecting miralayer interactions, it takes the analytical form 



T{q,iuJn,iu!m) = T^^tr {Q{iek)e"^''G{iek + iuj,n)^rg{iek + iuJn)e "''' + G{iek)e '^"Giiek - ii^n)^G{iek - icJ„i)e*'*'"} 
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FIG. 1: The diagrams contributing to the drag conductivity to leading order in the interlayer interaction f/(q, uj) (wavy hues). 
The fuU hues represent the electron Green function. The external vertices labelled by the velocity operator Vi are vector 
(current) vertices while the internal vertices are scalar (density) vertices. 
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FIG. 2: Diagrams defining the triangle vertex T{q,iJi,iJ2) 



where G denotes the Green function (for a particular reaUzation of the disorder potential) , e^ is a fermionic Matsubara 
frequency, and v represents the velocity operator. The vertex T should be averaged over reahzations of disorder, as 
vifill be discussed in Sec. IIIBI 

Summing over the Matsubara frequency w„, performing the analytical continuation to a real frequency fi, and 
finally taking the limit fJ ^ yields for the dc drag conductivitjii^*i& 
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167TTS ^ J_^ sm\?{uj/2T) 



rf^(q,w + iO,w-iO)rf (q,u;-iO,w + iO)|C/(q,w)|^ 



(5) 



In the sequel, we will use a short-hand notation, r(q, w) = r(q, w + jO,w — lO). Note that the Onsager relation 
ctIJ{B) = ajl{-B) implies, in combination with ((SJl, that r(q, w - iO, w + iO; B) = r(q, to + iQ,u) - iO; -B). 
The experimentally measured drag resistivity can be expressed via the drag conductivity as 
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where p\i^' are the resistivities of the layers. Note that the minus sign corresponding to the standard tensor inversion 
is absent in this expression, according to the conventional definition of the drag resistivity. This definition yields a 
positive transresistivity in the absence of a magnetic field. 

The triangle vertex T{q,uj) is obtained by analytic continuation of @, see Appendix IXI for detail. The result has 
the form T = F^"^ + F'^''^ with the two contributions 
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ctr {^rg- (e + u:)e'^'[g- [e) - 0+ {e)]e-''^'g+ {e + c^)} + (o., q ^ -uj, -q). 



(8) 



Here, G^{€) denotes the advanced/retarded Green function and fj, is the chemical potential. Note that at zero magnetic 
field only F^**) survives,^- whereas F^^^ containing products of three advanced or three retarded Green functions is 



zero. By contrast, in strong B both r*^°) and F^''^ should be retained. Most importantly, we will show below that in 
the ballistic limit there is a cancellation between F^"^ and F*^**) in the leading order. 
For small uj, the expressions for F(q, w) simplify to 

F('^)(q,c.) = ^M{vg+{e)e'^'^g+{e)e-^^'^g+{e) - {0+ ^ Q-)} (9) 

FW(q,'^) = -tr{vg-(e)e^i"-[g-(e)-g+(e)]e-"i'-g+(e)}. (10) 



For well-separated LLs, this approximation holds as long as to is small compared to the width A of the LL. It is also 
useful to note that T^''\q,uj) can be expressed as 



F('^)(q,c^) = -VqImtr{e''J'-^+(e)e-^'J'-5+(e)}, (11) 

which shows that F("^(q, a;) gives only a longitudinal contribution (parallel to q) to F(q, w). 

B. Impurity diagram technique in high Landau levels — SCBA 

In this subsection, we discuss the averaging over the random potential of impurities. We assume white-noise 
disorder, characterized by zero mean, (C/(r)) = 0, and by the correlation function 

(t/(r)W> = ^^(r-0, 

where j/q = m/27r denotes the zero-B density of states per spin and tq the zero-B elastic scattering time. We perform 
the averaging in the self-consistent Born approximation (SCBA) . This approximation, which neglects diagrams with 
crossing impurity lines, can be shown to give the leading contribution when the Fermi energy Ep is in a high LL 
with LL index N ^ li^ Strictly speaking, the disorder potential in the experimental samples is expected to be 
correlated on the scale of the distance of the two-dimensional electron layer from the donor layer. However, we find 
that the experimental observation can already be understood when considering white-noise disorder and that a finite 
correlation length of disorder does not qualitatively change our conclusions .-iS, 

Within the SCBA for well-separated Landau levels, ^^ the impurity average of the Green function, denoted by G^(e), 
is diagonal in the LL basis \nk) in the Landau gauge and takes the expression 

with the LL energies £"„ — ujc{n + 1/2). For energies e within a Landau level, the self-energy is given by 

S^(6) = lie -E^± z[A2 - (e - Er.^/'}. (13) 

Here, the LL index n is chosen such that \e — i?„| < A. The LL broadening A can be expressed in terms of the 
zero-field scattering time tq as 

A^ = 2uJe/TTT0. (14) 

The density of states is 

iy{e) - 1/^2^2^V(6) = i.oTo[A2 - (e - i?„)2]i/2 (15) 

with the elastic scattering time 

T{e)^[A'-ie-E,,r]-'/\ (16) 

Here, i = (l/eS)^/^ denotes the magnetic length. 

In principle, disorder leads to vertex corrections of both the vector and the scalar vertices of the triangle diagram 
F. However, for white-noise disorder there are no vertex corrections of the vector vertex. The vertex corrections of 
the scalar vertices generally involve impurity ladders (cf. Fig. O and turn out to be independent of the LL indices n 
and n'fSi 

l!::^,,,{e + u;,e;ci) = r^i^,^) {nk\e^'''Wk'). (17) 
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FIG. 3: Diagrammatic representation of the equation for the vertex corrections 7^^^ ^,y (e + u, e; q) (full triangle at vertex) of 
the scalar (density) vertices in the SCBA. Dashed lines represent impurity scattering. We also indicate that for well-separated 
LLs, the internal Green functions in the right-most diagram should be evaluated in the valence LL A'' which can differ from the 
LL labels n, n of the external Green functions. 

Here, the indices ^^v — ± indicate the type of Green functions involved in the vertex. In the limit of well-separated 
Landau levels, one finds the explicit expressions for the vertex corrections at w = 



7++(q,c.) = 



1 



l-j2(g^^)[^/2£-]2 
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(,q,^j 



1 - J^{qRc) 



(18) 
(19) 



where Jn{z) denotes the Bessel functions. The derivation of these expressions is reviewed in Appendix IbI 

For later reference, we also collect relevant matrix elements between LL eigenstates \nk) in the Landau gauge 
A = B{0,x). The vector vertex involves the matrix elements 



{nk\v^\n'k') = Skk' 
{nk\vy\n'k') = 5kk' 



meV2 

1 



m 



eV2 



{y/nSn^,i' + l - Vn + lSn,n'-l} 
{VnSn,n' + l + Vn+ lSn,n'-l} ■ 



(20) 
(21) 



In the limit of high Landau levels, n ^ N ^ 1, one can use quasiclassical approximations for these matrix elements, 
namely (nfc|ua;|n ± Ik') ~ ^iSkk'Vp/'i and {nk\vy\n ± lA:') ~ Skk'Vp/'^, with the Fermi velocity vp. 
The scalar vertex involves the matrix element 



{nk\e'^'^\n'k') = S.^^k^k' 



■ exp-T?^ 

n\ 4 



qa:ik + k')t][{qy + iq^)£'_ 



n — n' jn — n /^2/j2 



L 



{q^t/l) {n > n'), (22) 



where L"^ is the associated Laguerre polynomial. The expression for the matrix element for n < n' can be obtained 
from (|22|l by complex conjugation with the replacement q — > —q, nk <-^ n'k' . Since the characteristic LL indices are 
large, n,n' ^ l,\n — n'\, and relevant momenta are small compared to the Fermi momentum, q <^ kp, Eq. H22fl can 
be simplified by using the quasiclassical approximation. 



{nk\e'^''\n'k')~Sg^^k-k'f' 



-i0g(n-n')g-J<?x(fc+fc')^^/2 J _ ,((7/j(™)) 



(23) 



where (pq is the polar angle of q, i?c = i\/2n -I- 1 is the cyclotron radius of the n-th LL, and m = (n -I- n')/2. For 
most of the calculations below, the dependence of the cyclotron radius on the LL index in the vicinity of the Fermi 
level will be immaterial, and we will drop the corresponding superscript and simply write Re- The n-dependence of 
Re will, however, be crucial for the evaluation of the contribution to the drag related to the curvature of the electron 
spectrum, see Sec. IIII C 31 In view of the rotational invariance of r(q, w), it is sufficient to calculate it for a certain 
direction of the wave vector q. Choosing q to point along the positive x-axis, we simplify (|23|l to the form 



{nk\e''i-\n'k') ~ 4fe'*"-"'e-'«"V„_„-(gi?c). 



(24) 
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FIG. 4: The diagrams contributing to the triangle vertex in SCBA to leading order in the limits of well-separated Landau 
levels A/oic and large A^. 



III. TRIANGLE VERTEX r(q,w) 
A. Leading order 



We now turn to an evaluation of the disorder-averaged triangle vertex r(q, w) for well separated Landau levels, 
A/o^c <C 1, in the limit in which the Fermi energy is in a high Landau level, N ^ \. The relevant diagrams are shown 
in Fig.^ We begin by considering the low-temperature limit, cl',T <^ A. 

In the limit under consideration, the calculation is simplified as follows. Very generally for white-noise disorder, 
there are no vertex corrections of the vector vertex. By contrast, vertex corrections of the scalar vertices have to be 
retained. To leading order in A/wc, two of the three Green functions in Eqs. lO and (|10() should be evaluated in 
the A^th Landau level in which the Fermi energy is situated. Since the velocity operator has matrix elements only 
between states in neighboring Landau levels, one of the Green functions adjacent to the vector vertex must be taken 
in Landau levels A^ ± 1. This is illustrated in Fig. 0] 

We first consider the contribution r(")(q, w). In this case, it is most convenient to start from the simplified 
expression in Eq. (|ll|l in which to leading order in A/wc, both remaining Green functions can be evaluated in Landau 
level A^. Using the identity Vq[Jo((i'i?c)]^ = —'2(lRcJo{QRc)Ji{qRc) (with q = q/q), one obtains 



r(")(q,c^) - -2q^ Jo{qRc)MqRc 



{[G^7 



[GnI-'' 



f}, 



(25) 



where 7'^" = 7^''(q, 0) and the factor 2 accounts for the spin degeneracy. The calculation for r('')(q, w) in Eq. H1U|) 
yields 



r(^)(q,c.) = 2ci^ Jo{qRc)JliqRc) ^ [G+ 7^ 
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Summing both contributions, one obtains 



MqRc)JliqRc) Rc [G+ (7++ - 7+")] Im [G+7 



(26) 



(27) 



For arbitrary T,u! < u!c, this contribution takes the formSS 

r(q,w) =^ " 



de 



tanh- 



2T 



Re [7"+(q, uj) - 7++(q, w)] Im7++(q, w). 
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e-/x 



2T 



(28) 



Since the interlayer interaction is suppressed at large momenta q by a factor e~'°, where a is the interlayer distance, 
the drag conductivity ((SJ is governed by momenta q < 1/a. Depending on the relation between a and the cyclotron 
radius Rc, one distinguishes between the diffusive (a ^ Rc) and the ballistic (a ^ Rc) regimes. While in the former 
case, only "diffusive" momenta {qRc ^ 1) are relevant, in the latter case both "ballistic" {qRc ^ 1) and diffusive 
momenta contribute to the drag conductivity ^. Experimentally, when the transresistivity is measured in moderately 
strong magnetic fields (i.e. in high Landau levels), the condition Rc > a is typically satisfied. For this reason, we 
mainly concentrate on the ballistic regime in this paper. In Sees. IIII Bl and llil CI we will calculate the triangle vertex 
r(q, w) in the diffusive and ballistic ranges of momenta, respectively. These results will be used in Sec. |3for the 
calculation of the drag resistivity. 



B. Diffusive momenta, qRc <C 1 

In the diffusive range of momenta, qRc <C 1, we can expand the Bessel functions in the expressions for the vertex 
corrections, Eqs. (|18|) and H19|) . Due to the singular behavior of the vertex correction 7"* at small momenta g, we 
have 7"' 3> 7^^, 7 , so that only the contribution proportional to 7"' should be retained in (|27|l . This yields 

More generally, at small momenta one should also take into account the frequency dependence of 7"' , which has the 
structure of a diffusion pole, 

7+~(q,cj) = — (30) 

T[e)[D{e)q'^ — icjj 

where D{e) ~ i?^/2T(e) is the (energy-dependent) diffusion constant in a strong magnetic field. Eq. if^ is then 
generalized to 

„ AbjqRl D{ti)q^ , , , , 

n^^-) = ~^^^.jU(i;^^A^-E^)- (31) 

This result can be recast in the form (rig is the electron concentration) 

er,(q,c^) = 2-^ • g, Imn(q,c^) ~ 2-p^ ■ q, Imn(q,c^), (32) 

which allows for a simple interpretation as a nonlinear susceptibility^^^ This rewriting of Eq. (|31|) uses the resulli^ 
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for the diagonal conductivity in SCBA [daxx/due ^ doxy/drie for separated LLs) and 



(33) 



n(q,c.)^2KM) n,^y . , (34) 

for the polarization operator in the diffusive limit. It is worth emphasizing that the diffusive result in Eq. H31() 
arises from F^''^ only, since the other contribution T^"'^ does not contain the vertex correction 7^ . Note that the 
authors of Ref. Il2 failed to obtain the leading diffusive contribution H31|l . (|32|l . because of an incorrect treatment of 
vertex corrections2i*2^. For the same reason, they missed the 0{l/qRc) contribution [Eq. (|36|l below] which becomes 
important in the ballistic regime, as we are going to discuss. 

C. Ballistic momenta, qRc 3> 1 

1. Cancellation of leading contribution and the 0{l/qRc) contribution from vertex corrections 

In the leading order in the ballistic range of momenta, qRc ^ 1, the vertex corrections in Eqs. (|18(l and (|19|l can 
be neglected, 

7++(q,c^)~7+-(q,c^)~l. (35) 

Inserting this into Eq. 1)27(1 for the triangle vertex r(q, w), we immediately see that the triangle vertex vanishes to this 
order. We emphasize that r'"^) and F'''' do not vanish separately but rather cancel each other in the leading order. 

Thus, to obtain a non-zero answer for F(q, lo) from Eq. H27|l . we need to consider the vertex corrections in Eqs. 118|) 
as (I19|l in next-to- leading order in qRc- In this way, one finds from Eq. (|27|l 



(\/\y 



n±l 




(\/\y 




nil 



FIG. 5: Diagram contributing to corrections of order A/cJc to tlie triangle vertex. Here botli Green function adjacent to the 
vector vertex should be evaluated in Landau levels different from A'^. 

Here, we introduced a superscript on r(q, cj) in order to distinguish this contribution from other contributions 
computed below. At finite T and lu (assuming T,lu < lOc), wc find 
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(37) 



The contribution H36|) to r(q, cj) has been obtained in leading order in the limit A/uJc <C 1 and q/kp ^ 1 and in 
next-to-leading order in qRc 3> 1 . Thus, we are also forced to consider separately next-to-leading order corrections in 
the parameters A/ujc <C 1 and q/kp ^ 1, with the other two parameters kept in leading order. 

Before we turn to these calculations, we briefly remark that the leading-order cancellation in the ballistic regime 
was missed in Ref. '14 since the contribution from T^"^ was overlooked. The results obtained there for the diffusive 
regime remain valid since in this case F^"^ is negligible compared to T'^''\ see Sec. 1111 J^ 



2. Contributions of order A/cuc 



In this section, we consider the first corrections to the leading order in A/cOc to the triangle vertex r(q, w), while 
working to leading order in the ballistic limit qRc ^ 1 for high Landau levels iV ^ 1. While such corrections are of 
higher order in the small parameter A/lOc, this smallness may be compensated by a large factor qRc since it turns 
out that in this case there is no cancellation between F*^"^ and F*^**^ . 

Corrections of order A/ujc arise from two sources: (i) The Green functions adjacent to the current vertex are both 
evaluated in Landau levels different from N. (Note that the Green function between the scalar vertices must still be 
evaluated in the iVth Landau level because G^ — G~ ~ A/w^ for n ^ N.) This contribution is depicted in Fig.[Sl (u) 
The diagrams in Fig. 0] can be evaluated more accurately, keeping corrections in A/uJc, which arise from keeping the 
self-energy parts of the Green functions of Landau levels A^ ± 1. Note that we may now neglect vertex corrections at 
the scalar vertices because we consider the leading order in qRc 3> 1. 

Details of this calculation are presented in Appendix[ni Here we only state the results. The contribution (i) vanishes 
for both F'^°) and F'^''^. The contribution (ii) turns out to still give a vanishing contribution to the longitudinal triangle 
vertex, due to the cancellation between F*^"^ and F^''^ described above. However, the transverse contribution to F^**) 
no longer vanishes when considering corrections in A/lUc- In this way, we obtain the contribution 
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(38) 



to the triangle vertex. 

As we will see below, the A/ujc contribution is of crucial importance for understanding the experimental findings. 
We mention that this term was lost in Ref. [13 (in addition to the l/qRc contribution missing there because of 
an inaccurate treatment of vertex corrections) in the course of the so-called "triangles-to-bubbles" transformation. 
Specifically, in Ref. [ij the self-energy in the Green functions connected by the current vertex was neglected compared 
to the cyclotron frequency, which obviously misses corrections of order of A/lOc. 

Eq. 1)38(1 is derived in the low-temperature limit, when T,uj ^ A. To analyze the temperature dependence of the 
drag, we will need the A/wc-contribution also at higher temperatures. We find for arbitrary relations between T and 



and T,uj < LJc that 










r(^/-=)(q,c^) = -qxz-^Ji(gi?,)Jo(gi?c) 




X / de tanh ^ - tanh 

J-oo \ 2T 2T J 




(-En ^ 
X Re 

A 


r {e + uj/2-ENf] 
[ A2 J 


1/2 

Re 


[l (- 


- a;/2 - En? 

A2 



-,1/2 



(39) 



5. The conventional eontribution of order q/kf 



In this section, we compute the contribution to T due to terms of order q/kp relative to the leading order. Such 
terms arise from a more accurate treatment of the matrix elements involved in the scalar vertices, for which we now 
use the more accurate expressions 



{N±l\e'^^\N){N\e-'^^\N) ^ i,h{qR,[l± ^])UqR,) 



(iV|e'i""|^>(^|e""'''|A^±l) 



lJ^{qR,[l±—])UqRc) 



together with 



J,{qR,[l ± ^]) c. Ji(<zi?e) ± ^MqRc)- 



(40) 
(41) 

(42) 



Thus, such terms give rise to a contribution of the order of q/kp relative to the naive leading order [which vanishes 
because of the cancellation between F^"' and r^''^]. 

We find that such corrections arise only for the contribution T^^' , yielding 



r(9/fe^)(q,tj) =qx z — 



8w A2 ~ (/i - En J ,2 



J^iqRc). 



TT- A4 

Similarly to Eq. H39|l . we generalize this 0{q/kp) contribution to the finite-T case. 



(43) 



^{q/kf 



'(q,^) 



q X z 



X Re 



4Jo'(9^c) 



2A2 



7r2A 



, / -, e + uj/2 — n , e — uj/2 — n 

de tanh - tanh - ' ^ 



2T 



{e + uj/2-ENf 



2 1 1/2 



A2 



Re 



(e - uj/2 - En) 



2T 

21 1/2 



A2 



This expression can also be rewritten as 



2a T. 



er(9/'=^)(q,w) =qx z ^^ Inin(q, tj) 



(44) 



(45) 



with the polarization operator n(q, lo) for the ballistic regime [cf. Eq. H5U|) below] and the Hall conductivity 



eue e2 ^^ A 
B 7r2 wc 



1- 



(a* - EnY 



3/2 



(46) 



in SCBA. It can be checked that Eq. H45|l is valid for arbitrary T, including T > loc- 

The q/kp contribution arises from taking into account the dependence of the cyclotron radius and hence the velocity 
on the Landau level number, which is a direct consequence of the curvature of the zero-B electron spectrum. It is thus 
natural that the obtained result (|45l) is a high-magnetic field analog of the conventional contribution to FM Only this 
contribution was retained in Refs. lZil3t while the other contributions related to the particle-hole asymmetry due to 
the LL quantization were lost there. 
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IV. SCREENED INTERLAYER INTERACTION 

In this section, we summarize the resuhs for the screened interlayer interactioni^ 

Vi2{q) 

Here, V^(q) = 27re^/(j denotes the bare miralayer interaction and V^iq) — V{q)e^'^°' is the bare mterlayer interaction, 
a denotes the distance between the layers. The polarization operator of layer I is denoted by 11; (q, a;). For q small 
compared to the Thomas-Fermi screening wave vectors kq,; = 47re^i'o,i (^ = Ij 2 labels the layer and I'o,; denotes the 
zero-field density of states per spin of layer I), this can be approximated as 

Ko,i'^o,2 smh(ga) ni(q,w) 112 (q,w) 
In the random-phase approximation, the polarization operator in a strong magnetic field has the form 

n(q,c.) = ^E^"-™(9^^)/ ^^fW{G+(e-t-c.)[G+(e)7++(q,c.)-G;;,(6)7+-(q,c.)] 

+ G;;(e-^)[G+(6)7+-(q,c.)-G„(6)7-(q,c.)]}, (49) 

where np(t) — 1/{1 + exp[(e — [i)T\\ = {1 — tanh[(e — /^)2T]}/2 is the Fermi distribution function and we have used 
the quasiclassical approximation for matrix elements H24|l . 

We turn now to a brief summary of results for n(q, cj) in various relevant domains of momenta and frequency. 
Some of these results can be found in Ref. [l^ we reproduce them here for the sake of completeness. The polarization 
operator in the diffusive regime of momenta (and at T ^ A) was already given in Eq. (|34|l . In the ballistic regime 
qRc ^ 1, the expression H49|l can be simplified by neglecting the scalar vertex corrections, 

1 t°° de 

n(q,c.) = —Y^jl_^^{qR^)j^^-nF{e) [G+(e + c) + G-(e - ^)] Im G+ (e), (50) 



For low temperature and frequency, w,r ^ A, the real part of the polarization operator H50() takes the form 



Rcn(g>i?-i,cj^O) = 2vo + 2vo^J^{qR,) 



^_{^,-E^^'"' 



A2 



(51) 



Here, the first termS^ arises from Landau levels with n ^ m, while the second term represents the contribution of 
the A^th LL {n = m = N). The intra-LL (second) term contains an additional energy factor |[1 — (/i — E^Y/A'^] 
compared to the case of diffusive momenta, which is due to the suppression of vertex corrections at high momenta. 
The imaginary part of the polarization operator for w, T ^ A has the form 



ImH(q,c^) = 2i.o^Jo('Z^c) 



^ {f^-E^y 



(52) 



A comparison with Eq. H51|) shows that ImH ^ ReH in this regime. 

It follows from Eq. H51|) that there is an additional wavevector scale qRc ~ uJc/A in the ballistic regime, where the 
behavior of ReH changes. Specifically, for q <C oJc/ARc the polarization operator (and hence screening) is due to the 
contribution of the A^-th Landau level, while at larger q it is due to Landau levels with n ^ N. Only in the latter 
case, we recover 

H(q,c^)~2i.o, (53) 

and thus the standard B = form of screening. When the temperature is large compared to the Landau level 
broadening, A ^ T -^ uic, Eqs. (|34|) . (|52|) . and the second term of (|51|l are effectively multiplied by factors ^ A/T 
due to thermal averaging. In this case, the real part of H(q, w) takes its zeio-B form under the weaker condition 
qRc S> LOc/T. This follows from the expression 



ReH(q,c.) =2.0 + 2.0^^ j'oilRc)Q{^) cosh"^ F^^ • (54) 
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Here we defined the function 



2 1/2 , _ sgn{z + 2x) _ /2 _ sgn(z-2x) _ ^ _ ^/^ 



1 



Q(x) = / dzz{\~zy\z--^^ ^Re[(z + 2x)'-l]^/^-^^- ^Re [(^ - 2x)" - 1]^/" ^ . (55) 



Tlie imaginary part for T <^lOc reads 



Im n(q,c.) = 2.of^Jo^(gi?c) ^7^ (^) cosh- (^i^ j , (56) 

where 7i(x) is a dimensionless function representing the overlap of two Landau bands, 

/oo 
dz JRe [1 - (z + x)2] '/'I JRe [l - (z - x)^] '^H . (57) 

Finally, in the high-T limit, T ^ Wc, the imaginary part of n(q, w) becomes independent of En — fj,, because of 
thermal averaging. 



Imn(q,w) ~ ^z^o— ^X! 



n.rn 



tanh '^"+'^~^' - tanh ^"'^ 
2T 2T 



Jn-,nilRc) n ( ^^^ 



= 2..i^X:4^„«,,H(^). ,58) 

Since 7i(|x| > 1) = 0, the imaginary part of n(q, w) as a function of tj at T ^ oj^ consists of a series of peaks 
(broadened by A) around multiples of the cyclotron frequency. 

Importantly, the imaginary part of the polarization operator is suppressed at high frequencies, lo 3> qvp- This 
follows from Eq. H58|) . since J^{qRc) is exponentially small when n ^ qRc- This is analogous to the zero-i? case, 
where 

lmU(c[,uj;B ^0) = 2i^o^^e(qvF-Lu) (59) 

qvF 

with 9{x) the step function, and can be traced back to the fact that at high frequencies the magnetic field becomes 
almost irrelevant, so that the polarization operator approaches its zero-i? formi^ 

V. DRAG RESISTIVITY 

In a strong magnetic field, uJcTq 3> 1, the intralayer Hall resistivity p^y dominates over the longitudinal resistivity 
Pxj;. Therefore, the drag resistivity is given by 

pS.^£^ <^yyP^'J- (60) 

Using Eq. (|SJ), we get the expression for the longitudinal component of the drag resistivity in a strong magnetic field, 

B B I f°° duj 



^'''' euei ene2 Stt J_^ 2rsinh^(w/2r) 

X y'^rW(q,c.,B)r(2)(q,c.,-B)|t/i2(q,c.)p. (61) 

The overall minus sign in Eq. H61|l is due to the relation p^y = —Pyx- It follows that for identical layers, the longitudinal 
(r ex q) component rn of the triangle vertex gives rise to negative drag, since r||(— B) = r||(— i3), while the transverse 
(r c>c z X q) component T ±_ yields positive drag, V^_{—B) = —r±{—B). 

Since the upper limit of the momentum integration in (|61|l is effectively set by the inverse interlayer distance, a~^, 
the behavior of the transresistivity will essentially dependent on the relation between Re and a. Below, we mainly 
concentrate on the ballistic regime 

Wc/A < Rc/a < NA/lu^, (62) 



which we consider as most relevant experimentally. In Sec. lVDl we will briefly consider other situations and discuss the 
evolution of the transresistivity with decreasing interlayer distance, from the diffusive {Rc/a <^ 1) to the ultra-ballistic 
{Rc/a ^ NA/lUc) regime. 
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A. Ballistic regime: Low temperatures (T <C A) 



In the low-temperature limit, the expressions derived for the triangle vertex r(q, cj) at w <C A are sufficient, because 
frequencies in Eq. I|51fl are restricted to uj ^ T ^ A. Let us analyze which of the contributions to the triangle vertex 
dominates, depending on the relation between q and 1/ Re- 
in the diffusive range of momenta, qRc «C 1, the leading contribution to the triangle vertex is given by Eq. H31|l : its 
magnitude can be estimated as 



r 



ujkp 
^qRc 



In the ballistic regime, qRc ^ 1, we have three competing contributions (see Sec. IIII (T| . 

p(A/^e) 



A2((ji?,)2' 






p(?/feF) 



A^Rc 



N 



r(i/9flc 



kp A 



(63) 

(64) 
(65) 
(66) 



Comparing these expressions, we ffiid that the first contribution, r'^^/'?^'^\ dominates for qR^ ^ ^c/A, the second 
one, T'^^/^''\ is dominant for lOc/ A <C qRc ^ NA/uJc, while the last contribution, r'^^/'"'^\ becomes the largest one 
for qRc ^ NA/iOc- This is valid provided that the Landau level index N is sufficiently large, N > (wc/A)^. We will 
assume below that this condition is fulfilled. 

Splitting the momentum integral in Eq. (|51|l into three parts, corresponding to regions of different behavior of the 
triangle vertex and polarization operator, we present the transresistivity in the following form, 



rxx 



kf 



87r3 V K^ne A2 ) ^ \ K^UcA"^ / 2 ^ -oo 



cLu- 



2Tsmh\u;/2T) 



/(w), 



I{uj) = hiuj) + Iu{lu) + Im{uj), 



(67) 



(68) 



where the subscript 1 — 1,2 in (. . .)j refers to the layer /, and the contributions /i. In, and /m in (|67|) are determined 
by the momentum domains qRc ^ 1, 1 ^ qRc ^ ^c/A, and qRc ^ ^c/A, respectively. The corresponding 
expressions are given in Appendix^ Estimating all three terms, we find [see Eq. (|D12(l ] that the leading contribution 
is given by the last term, 



I{uj) ~ /iii(cj) = 



1 



27r3a2i?2 



In 



RcA 



16 



{fi-EN)[A''-{fi-ENf 



16 



ifi - En)[a^ - {^i - ENf] 



1 



^i?2 



OJc 



RcA 



auJc 



Therefore, for T <C A we get for identical layers 

32 1 



rxx 



T\ fRcA 
3TT^e^ (kpafiKoRcy V A ' '' 



^i- En 



2 r 



1- 



it^-ENf 



(69) 



(70) 



Thus at low temperature T ^ A, the drag resistivity scales with the magnetic field and temperature as 

pZ^T^B\n{B,/B), (71) 

where _B* ~ {mc/e){vp/a'^TQ)^'^ sets the upper boundary for the considered ballistic regime on the magnetic field 



axis. 



If Re differs slightly between the two layers (i.e., the concentrations are slightly different) so that 6 Re/ a <C 1, the 
above calculation fully applies, with the only change in Eq. H7Q(I 



^i- En 



2 r 



1- 



{p - En)' 
A2 



N 



A 
H- En 



1 



1- 



{^i - EnY 



(Ai - EnY 



(72) 
(73) 
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FIG. 6: Low-temperature drag for identical layers. Left panel: temperature dependences of the C'(A/a;c)-term in p^^{T) for 
different values of the filling factor of the highest LL, i/jv — 0.5, 0.2, 0.8 (from top to bottom). Right panel: dependence of the 
0{A/ujc)-teim on the filling factor vn for different values of temperature, T/A = 0.1, 0.5, 1 (from top to bottom). 



This yields an oscillatory sign of the drag. For identical layers the drag is positive, at variance with Ref. Il4i This is 
because the leading term here originates from the component r_L of the triangle vertex transverse to the wave vector 
q (i.e. directed along z x q). For a more detailed discussion of the sign of drag in different regimes, see Sec. IV Dl 
If a ^ SRc ^ RcA/ojc, the calculation still applies, but the argument of the logarithm changes, 



a LUc 



In 



Re A 
SRc u>c 



(74) 



B. Ballistic regime: Arbitrary T/A 



Having identified the leading contribution (coming from A/wc-term) to drag for temperatures small compared to 
the LL width A, we generalize the obtained result to the case of larger T (and correspondingly lo). As discussed in 
Appendix^ the only difference in the momentum integral in the A/wc-term is the replacement A -^ T under the 
argument of logarithm. Using Eq. (|39|l and assuming that the difference in Re between the two layers is not too large, 
dRc <C a, we express the C'(A/a;c)-contribution to the transresistivity as 



D\i^/'^-=) 



(Pf.) 



1 



{kFa)'^{KoRc 
duj 



In 



i?cmax[A, T] 



OLOc 



, 2Tsmh^Lo/2T) 
where !F{u!, jjl, T) is a dimensionless function of cj/A, {ji — En)/ A, and T/A 



[T{uj, ^l,T)]l[T{uJ, ^l,T)]2, 



(75) 



T{uj,fi,T) 



de f , e + Lu/2 — u , e — Lu/2—n 

— tanh - tanh 

A V 2T 2T 



E 



N 



{e + uj/2-ENf 
A2 



1/2 r 



(e - Lu/2 - E 



N 



A2 



1/2 



(76) 



For arbitrary T/A ^ 1 this can only be calculated numerically. In Fig. IHlwe present the results for the temperature 
dependence of the C'(A/ci;c)-contribution to drag as well as for its dependence on the filling factor v^ of the highest 
(partially filled, < i^at < 2) LL. It is worth mentioning that when temperature is varied at fixed filling factor (as in 
typical experiments), the chemical potential is varying as well, /i = ^{vn,T), which is taken into account in Fig. 

Consider the regime of temperatures large compared to the LL width, T ^ A. In this situation the LLs will be 
broadened by the temperature, so that typically /i — E^ will be of order T and thus much larger than A. Expanding 
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the tanh-terms in ifT^ in uj < 2A ^ T and e — En < A ^ T, we arrive at 
.^.(AM) 4 1 /i?^\ 

X fsinh^^^^cosh-^^^^^) fsinh ""^ ' ^ cosh'^ ""^ " ^ 
V 2T 2T J ^\ 2T 2T 

A \ 3 /.OO 7 

tJ y^^[-^(^/2A)]i[^(c.'/2A)]2, (77) 

where !F{x) is a dimensionless function similar to Eq. H57|) . It also describes the overlap of two shifted Landau levels, 
but has an additional factor [(e — _Ejv)/A]^ arising from the particle-hole asymmetry due to LL quantization, 

T{x) = f dz z^ JRe [l-{z + xf] '^'l JRc [l - (z - xf] '^'j (78) 

J —OQ 

The contribution (|77|1 scales as 

pf,(xT^3BV2in(r/i32). (79) 

Since the above C'(A/a;c)-term falls off quickly at T ^ A, we should analyze the contributions of the other terms. 
Let us first calculate the "conventional" term 0{q/kF), substituting Eq. H45(l in Eq. H61|l . Remarkably, the strong-B 
expression for the g/fc^^-contribution to drag resistivity reduces to the standard zero-B formiiiiSii&, 

/ D\(i/kF) 1 f°° du 

\Pxx) 



47re^neine2 J-oa 2Tsinh^(Lj/2T) 
X J ^q^[lmll{q,u)], [Imn(q,c.)]2 \Ui2{n,uj)\' . (80) 

Here all the information about the magnetic field is encoded in n(q, lj), Eq. (|49|l . 
For T 3> A, expanding the tanh-terms in Im n(q, to) just as before we find 

/ D\ii/kF) 3C(3) 1 /cjc\2 A 

[Pxx) - 2Tr^e^{kFa)^{KoRc)^\A) T 

X (cosh-^^^)^ (cosh-^^) J^^[H(W2A)],[7^(./2A)]„ (81) 

where C(x) is the Riemann zeta-function [C(3) ~ 1.202] and 7i(a;) is defined in Eq. H57(l . This contribution scales as 

The slower fall-off of the C'(q/fci?)-contribution (|81|) as compared to (|77|) can be traced back to the different nature 
of the particle-hole asymmetry underlying these two contributions. Specifically, the C'(A/a;c)-term (|81|) is governed 
by the particle-hole asymmetry due to the LL quantization. This is reflected by the factor e — E^ in (|76|) which 
after thermal averaging, translates into a factor in Eq. (|77|l which is asymmetric in En — /i. On the other hand, the 
"conventional" q/kp contribution is due to the curvature of zero-B spectrum and therefore is symmetric in e — En 
(and in En — A* after thermal averaging). In both cases the fall-off of drag at T ^ A is due to the absence of electronic 
states outside the Landau band (for |e — En\ > A). However, the thermal averaging of the odd function of e — En 
yields an additional factor A/T for each A/cJc-triangle vertex, at variance with the case of an even function of e — En 
determining C'(g/fcf)-contribution. 

Finally, we evaluate the contribution of C'(l/gi?c)-term. On the one hand, the thermal averaging suppresses each 
r(i/gfic) vertex by the factor (A/T)^, similarly to the C'(A/a;c)-term. This is again because of the particle-hole 
asymmetry due to the LL quantization. On the other hand, the peculiarity of the finite-T screening gives rise to a 
factor (r/A)2 in the momentum integral involving the C'(l/(7i?c)-term, see Appendix^ The remaining frequency 
integral yields the factor T/A, since the allowed frequencies are restricted by |ijj| < 2A <C T. As a result, the 
contribution of this term to the drag resistivity is inversely proportional to temperature for T 3> A, similarly to the 
conventional q/fci^-contribution. For simplicity we restrict ourselves to the case of identical layers, where we get 

.o.(i/,fl.) c(i/^^-) 1 A . ^^^ EN-fi ^En-JJ- .„„. 
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where c^^/i^") is a constant of order unity, 

^{i/qR.) ^ _ / dx[V{x)fW{x), (84) 

I" J~oo 

with the fmiction W{x) defined in Eq. (|D20(I and 

dzz{z + x) Re[l - (z + :e)2]1/2 J Ke[l - (z - a;)^]!/^ I . (85) 

-00 

This contribution scales as 

pf, (X -T-ii?^/^ (86) 

We thus conclude that the 0(l/(7-Rc)-contribution wins over the 0(A/a;c)-contribution for 

V a^c ) 
Comparing H81|) and (|83|l . we have 

io(i/9i?c)i (fc;^^"^/^^ Ittva j « ^' ^^^^ 

as foUows from Eq. H62|l . Therefore the C'(l/g_Rc)-contribution dominates the drag resistivity in the intermediate 
range of temperature. This contribution oscillates with changing the filling factor of the two layers; however, it is 
negative for matching densities, unlike the 0(A/cjc)-contribution. 

For higher temperatures, T > ojc, the terms related to the LL particle-hole asymmetry fall off rapidly due to the 
thermal averaging involving many LLs and thus the q/kF-term (i.e. the "conventional" contribution to the drag 
resistivity) soon becomes dominant. The drag resistivity is then always positive, independently of the difference in 
filling factors of the two layers. It monotonously increases with increasing T and takes the form 

^- - 'Mjk^ajh^^^i^S £^[H(./2A)].[H(./2A)]. 

1 ^^^''^^ocr2si/2. 



e2(fcj.a)2(Koa)2 \Ef J A 

This is almost the same result that is found in zero magnetic fieldii^i^ti^; the only difference is an extra factor 
^ Wc/A oc i?^/2 The reason for the emergence of the zero-B result is physically very transparent. Characteristic 
frequencies to ^ T ^ coc set a characteristic time scale T~^, which is much smaller than the time of the cyclotron 
revolution. At such times the electron motion is essentially unaffected by the magnetic field. The magnetic field 
enters, however, through the density of states v inside the LL, which determines the characteristic magnitude of 
Imll, and thus of T, see Eq. (|45|l . The w-integration in (|61|l thus results in an effective averaging of ly'^, yielding the 
factor LOc/A. It is worth mentioning that, for the same reason, the longitudinal resistivity pxx of a single layer is also 
enhanced by such a factor in the regime T ::^ cjc ^ A as compared to its zeio-B value, see e.g. Ref. ^^ 

For still higher temperature, T ^ vp/a, the quadratic-in-T dependence of the drag resistivity crosses over into 
the linear-in-T drag. This occurs because of the suppression of the imaginary part of the polarization operator 
[determining the q/fci^-triangle, Eq. (|45|l ] at lu ^ qvp, see Eq. (|59f) . As a result, the domain of w-integration is 
effectively restricted by w ^ vp/a (since q ;$ l/o), yielding the replacement T^ -^ Tvp/a as compared to the case of 
Wc ^ 7" < vp/a, 

pf. « ri?V2. 

Before closing this subsection, it is worth mentioning that in the above consideration we have neglected the contri- 
bution of magnetoplasmons to the drag (see Ref. 13 for details). While this contribution may become important for 
very high temperatures, T ^ ujc, it is negligibly small at relatively low T '^ A, which is the range of our main interest 
in the present paper. 
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FIG. 7: Schematic temperature dependence of drag in the ballistic regime for matched and mismatched densities. In the latter 
case the mismatch is chosen such that the drag is negative at low T (see text). Scaling of p^^ with temperature in different 
regions is indicated: I - Eq. ||nj, H - Eq. JZHJ, HI - Eq. (jHSJ, and IV -Eg. ll^ . 

C. Comparison with Experiment 

In this subsection we compare the results for the drag in the ballistic regime obtained above with experimental 
findings. We have found a sequence of different regimes of the temperature behavior of p^^, see Eqs. (|7T]| . (17^ . (jHSJl, 
(|88ll . All these results are schematically summarized in Fig. [T] The upper curves there depicts p^^ {T) for equal 
densities, whereas the lower curve corresponds to a mismatch in densities chosen in such a way that the Fermi energy 
is located in the upper half of the Landau band in one layer, and in the lower half in the other layer. As has been 
already emphasized, the drag at low temperatures is positive for matched and negative for mismatched densities. This 
sign of the oscillatory drag can be traced back to the fact that the dominant contribution to the triangle vertex is 
given by r^ ''^'=', which is transverse with respect to the momentum q. 

We now compare these results with a most recent and detailed study by Muraki et al}'^ of the Coulomb drag in 
the regime of high Landau levels. A comparison of our Fig. [T] with Fig. 3 of Ref. llOl reveals a remarkable agreement 
between the experimental findings and our theoretical results. In both the theory and the experiment, (i) p^^{T) 
shows a sharp peak at low temperatures; (ii) the sign of the drag in this temperature range oscillates as a function 
of the filling factor of one layer (at fixed filling factor of the other layer); (iii) the low-T drag is positive for equal 
filling factors and negative when the Fermi energy in one layer is in the upper half and in the other layer in the lower 
half of the Landau band; (iv) the high-T drag is always positive, independently of the difference in filling factors of 
two layers and increases monotonically with increasing T. Furthermore, it was observed by Muraki et al (see Fig. 2 
of Ref. llOl) that in the low-temperature regime of initial increase of p^^, as well as in the high-temperature regime of 
"normal" drag, the drag resistivity can be described by an empirical scaling law. 



Pxx oc 



B 



fiT/B). 



(89) 



Our results for the low-temperature, H71(l . and high-temperature, H88|l . increase of p^^ are in a nice correspondence 
with this prediction, with f{x) ~ x^. 

The magnitude of the low-temperature peak in the drag resistivity that follows from our theory also agrees with 
the experiment. Specifically, estimating Eq. H7U|I at T = 0.25A and [{p — £'jv)/A]^ = 1/2 by making use of typical 



Ko 



10« 



lO^^m, Re 



10 



"^m, we find p^^ 



10, in good agreement with 



experimental parameters, 
the result of Ref. llOi 

There is however a difference between our result (|71|) for the low-temperature scaling of drag and the interpretation 
of low-T data in Ref. 10. Specifically, while we find T^ scaling in this regime, Muraki et al. fit the data to an exponential 
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FIG. 8: Schematic temperature dependence of low-temperature drag in different regimes: a) diffusive, Rc/a <C 1; b) weakly 
ballistic, 1 < Rc/a < ujc/A; c) ballistic, ujc/A < Rc/a < NA/ujc; d) ultra-ballistic, NA/ujc < Rc/a. 



(activation-type) dependence, arguing that localized states are responsible for the low-temperature "anomalous peak" 
in Px^{T). We do not expect, however, that localization plays an important role in the regime of high Landau levels 
at realistic temperatures. Indeed, as is seen from Fig.l of Ref.llOl the resistivity for filling factors ly ^ 10 has a shape 
as predicted by SCBA, without developed Hall plateaus. Also, the fit of the low-T behavior of p^^ to the activated 
over a single decade is not unambiguous; the same data could be quite well fitted to the T^ power law. In other words, 
we believe that our theory based on SCBA and not including quantum localization effects is sufficient to explain the 
most salient experimental observations of Ref. llOt the "anomalous" drag with oscillatory sign at low temperatures 
and the "normal" positive drag at high T. 



D. Evolution of p^^{T) with varying interlayer distance: From the diffusive to the ultra-ballistic limit. 



As discussed in the beginning of Sec. lV Al the form of the drag resistivity p^^{T) depends on the value of the ratio 
Rc/a. In the above we have concentrated on the regime uJc/A <C Rc/a <C NA/lUc, which can be termed "ballistic" 
and which we believe to be most relevant to a typical experiment. In this subsection we briefly describe the results 
obtained for other regimes. Specifically, with increasing Rc/a we identify the following four regimes: i) diffusive, 
Rc/a <^ 1, ii) weakly baUistic, 1 <C Rc/a <C cOc/A, iii) ballistic, uJc/A <^ Rc/a <^ NA/uJc, and iv) ultra-baUistic, 
NA/ujc ^ Rc/a. In all regimes, the temperature-dependence of the drag resistivity is non-monotonous: the absolute 
value of Pxxi"^) shows a peak around T ^ A and increases again at T ^ Wc- However, the T— and B— dependences 
oi p^^, as well as the sign of the low-temperature peak (the high-temperature drag is always positive), are specific for 
each particular regime, as illustrated in Fig.|Hland summarized below. 

Diffusive regime, Rc/a <^ 1. In the diffusive regime, the drag at not too high temperatures, T <^ ujc, is governed by 
the diffusive rectification, Eqs. (|31|l and H32|l . As a result, the sign of the drag at T ~ A oscillates but is opposite to 
what we found above for the ballistic regime: the drag is negative for equal densities>i^ At the "slopes" of the peak, 
p^^ scales with T and B in the following way 



Pfx« 



- r2 ln(rB3/2), 

- r-iB3/2 InB, 



r< A, 

T> A, 



(90) 



where the sign corresponds to the case of matching densities. 

Weakly ballistic regime, 1 <C Rc/a <^ lOc/A. This regime is qualitatively similar to the diffusive regime. The peak 
at T '--^ A is governed now by the 0{l/qRc)-teTia in the triangle vertex, resulting in 



r< A, 

A < T < T„ = LUc{a/Rc) 
T:^uc{a/Rc), 

The sign of the peak oscillates just like in the diffusive regime. 




(91) 



18 



p-h symmetry p-h asymmetry 

curvature LL - quantization 




I^ 



p-h symmetry 



o 

h O 
O- 



Je= -Ji 



p-h asymmetry 

curvature DoS V(E) 

V V 



O- 

o- 



o— 



o- 
o- 
o- 
o- 



Je^ -Jh 



FIG. 9: Schematic illustration of different sources of particle-hole asymmetry: curvature of zero-B spectrum E{k) vs LL- 
quantization of the density of states (DoS) v{E). In the particle-hole (p-h) symmetric case, the electronic and hole contributions 
to the current induced in the passive layer {j^ and jh, respectively) compensate each other. When the p-h asymmetry is generated 
by a finite curvature, the velocities of electrons and holes (shown by arrows in the right panel) are different, which destroys 
the compensation. This is the "conventional" mechanism of the drag. When the DoS depends on energy (in the present case 
because of the LL-quantization), an "anomalous" drag arises due to the difference in numbers of occupied electronic and hole 
states. 



Ballistic regime, Wc/A <C Rc/a <^ NA/loc. This is the regime we have studied in the main part of the paper. For 
the reader's convenience, we repeat the results here. The peak is governed by the C'(A/a;c)-contribution, its sign 
oscillates and is positive for matching densities, 



n 



T2 B \n{B,jB), 
oc { T-3 5^/2 \n{B,lB), 



r< A, 

A<T<T,, 



(92) 



Ultra-ballistic regime, NA/lUc <^ Rc/a- The drag for all temperatures is determined by the conventional 0{q/kp)- 
contribution and is always positive. 
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(93) 



At high temperature, T ^ tOc, the drag is governed by the conventional contribution (and is therefore positive) in 
all the regimes. It is linear in T in the diffusive regime {p^^ ex TB"-^/^). In all the ballistic regimes the drag resistivity 
scales as pf^ ex T'^B^/'^ for ujc^iT <^ vp/a and pf^ oc TB^/"^ for T > vpla. (As mentioned in the end of Sec. lVBI 
we do not consider the magnetoplasmon contributiorJ^ here.) 



VI. SUMMARY 



In this paper, we have developed a systematic diagrammatic theory of the Coulomb drag in moderately strong 
magnetic fields, when the Landau bands are already separated but the Landau level index is still large. Using the 
self-consistent Born approximation, we performed a thorough analysis of all relevant contributions and, on this basis, 
analyzed the temperature dependence of the drag resistivity. Depending on the relation between the cyclotron radius 
Re and the interlayer distance a we distinguish several regimes. We concentrated on the experimentally most relevant 
ballistic regime. In this case the theoretical analysis requires special care, in view of a cancellation between leading- 
order contributions to the triangle vertex T. We also briefly considered the evolution of the drag resistivity in the 
whole range of Rc/a, from the diffusive to the ultra-ballistic regime. 

We have shown that Coulomb drag in strong magnetic fields is an interplay of two contributions, as illustrated in 
Fig. 1^1 At high temperatures, the leading contribution is due to breaking of particle-hole symmetry by the curvature 
of the zeio-B electron spectrum. This "normal" contribution to the drag is always positive and increases in a broad 
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temperature range as T^. At low temperatures, we find that a second, "anomalous", contribution dominates, which 
arises from the breaking of particle-hole symmetry by the energy dependence of the density of states related to Landau 
quantization. This contribution is sharply peaked at a temperature T ~ A (where A is the Landau level width) and 
has an oscillatory sign depending on the density mismatch between the two layers. In particular, we find that in the 
ballistic regime the sign is positive for equal densities, in contrast to the negative sign in the diffusive regime found 
in Ref. Hi 

Our results for the temperature dependence and sign of the drag resistivity p^x{T) in the ballistic regime are 
illustrated in Fig.EI These results are in good agreement with recent experimental findings^", and thus explain the 
remarkable features of Coulomb drag in high Landau levels observed experimentally. 

Finally, we discuss some prospects for future research. First, our theory can be generalized to phonon drag, which 
is expected to dominate over Coulomb drag at larger separations between the layers. Second, it will be interesting to 
consider the magnetic field and temperature dependence of the drag around filling factor v = 1/2, where transport 
is due to composite fermions moving in a reduced magnetic field^ Third, one can study the effects of quantum 
localization, as well as criticality in the center of the Landau band^^ which should become important in lower Landau 
levels or for very low temperatures. Finally, it should be possible to reproduce our results within the framework 
of a quantum kinetic equation [cf. Ref. 24] . This would also allow one to generalize the theory of magnetodrag to 
non-equilibrium setups (strong bias, microwave, etc.), as well as to other observables (e.g., the thermopower) related 
to particle-hole asymmetry. 
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APPENDIX A: ANALYTICAL CONTINUATION 

In this Appendix we perform the analytical continuation of the Matsubara expressions for the drag conductivity and 
the triangle vertex T. To calculate the Matsubara sum over uj„ = 2TmT in Eq. (|3J), the standard contour integration 
in the complex to plane is done, 

Tj2fi^^n) = ^J dlO f{L,) COth^. (Al) 

The integrand has branch cuts at Im w = and Im u) = —^k, where fife represents the external frequency. The 
integration contour Cf, thus contains three parts, see Fig. ^] Deforming the contour as shown in Fig. 1101 we get four 
terms corresponding to four lines (above and below of both the branch cuts) forming the new contour. 



g2 roo ^ 



rf ^ (q, u; + ink,uj + iO)rf '(q, LJ + iO,uj + iOfc)C/(q, lu + if7fc)C/(q, lj + iO) 
~ rf (q, UJ + i^k, UJ - iO)rf ^(q, uj~iQ,LO + iVLk)U{(i, u + if^fc)C/(q, lo - iO) 
-K rf ^ (q, uj + iQ,uj- i^k)^f (q, uj-ink,uj + iO)C/(q, uj + iO)C/(q, uj - iVtk) 
- rf^(q,t^-iO,cj-if7fe)rj^^(q,a;-irjfe,cj-iO)C/(q,cj-iO)C/(q,cj-ir2fc) . (A2) 

In the third and fourth terms we have used coth(z-|-iJ7fc/2T) = cothz. The contributions of points w = and w = —iVlk 
cancel the integral over the small circles around these points, so that the integrals above should be understood in the 
principal value sense. 
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FIG. 10: Contours for the tj-integration. 
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FIG. 11: Contours for the e-integration. 



We now perform the analytical continuation iflk ^ il + iO and take the limit fJ ^ 0. As shown in Ref. iWL the first 
and the last terms coming from outer sides of branch cuts vanish in the limit f7 ^ 0. This yields 



e 



E 



cLu coth 



U! d 

2T 'K) 



rf ^ (q, uj + i(),uj- iO)rf' (q, uj-iQ,uj + iO)C/(q, uj + iO)[/(q, uj - iO) 



Using 



d ^ UJ 1 

T^coth^— = 7^ , 

duj 2T 2Tsmh^{uj/2T) 



(A3) 



(A4) 



we arrive at Eq. (5). 

The next step is the analytical continuation of the triangle vertex. The summation over the fermionic Matsubara 
energies e^ = 27r(fc + 1/2T in Eq. Q) is performed using the integral 



7^E/(*^^) = 7^ / rfe/(e)tanh 



€ 

2T' 



(A5) 



along the contour Cf shown in Fig. 1111 Since the triangle vertex depends on two frequencies iujm and iaj„, the 
integrand now has three branch cuts in the complex plane of e, namely at Im e = 0, Im e = —uJrm and Im e ~ ~0Jn- 
Similarly to Cf,, the contour C/ can be deformed into a set of six lines going on both sides of each of the branch cuts 
(see Fig. ^fl j yielding 



f°° de , e 
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X tr{v [g+{e)^g-ie)y^''gie-iUn)e-'''''g{e + iUrn-iLOn) 

- V g(e + iuj,,)e'^''[g+{e) - g- {e)]e-'^''gie + luo^) 

+ V ^(e - luj^ + iuj,,)e'^'g{e - iuj^)e-'^'^[g+{e) - g- [e)]] 

+ (w„ -^ -ujm, q -^ -q)- (A6) 



In this formula ^ (e) — tj(e±i0) and we have used tanh(z — iWrn/ST) = tanh(2; — «a;„/2T) = tanhz. The equation (|A6p 
is vaUd irrespective of the relation between uj^m w„, and 0. Performing the analytical continuation to real frequencies 
iojrn -^ uji + iQ and iuon — > ^2 — iO (and shifting the integration variables e — > e + Cij2 a-nd e — s- e + o^i in the first and 
third terms, respectively) we obtain 



r(q,., + ^o,.2-^o) = y_- 



X tr jtanh ^-^ v [g+{e + UJ2) - ^"(e + u^2W'g+{e)e-''''g+{e + c^i) 

- iB.n\i^wg-{e + u:2)e'^''[g+{e)-g-{e)]e-'^'g+{e + u:^) 

+ tanh ii^ V g-(e + c^2)e*'i""e-(e)e-"i'-[g+(e + Wi) - g- {e + o-i)]} 

+ (cj,q^-cj,-q). (A7) 

Setting uji ~ uj2 and collecting the contributions containing only retarded (from the first term) and only advanced 
(from the third term) Green functions, we arrive [up to a redefinition of zero of fermionic energies, which are counted 
from the chemical potential in Eq. (|A7(l ] at Eq. Q for T^°-\ The remaining terms constitute the expression for 

APPENDIX B: VERTEX CORRECTIONS IN SCBA 

In this appendix, we review vertex corrections in SCBA. We start by noting that in real space, the impurity-averaged 
electron Green function in SCBA can be written as 

G(r, r'; E) = e^'^^""''-') J] C„(r - r')Gn{E) (Bl) 

n 

with 

G.(.O^^e-."'.V-i..(<L^)^ (B2) 

The gauge-dependent phase (p{r, r') satisfies ip{r, r') = —ip{r', r). This can be used to express the vertex correction in 
real space as (cf. Fig. EJ 

7'^''(q, uj; r) = e'^^" + -^ — / rfr'7^''(q, w; r')G^(r, r'; E + w)G'^(r', r; E). (B3) 

For well-separated Landau levels, the valence LL with LL index N gives the dominant contribution so that 

r^in, Lo; r) =. e^^^ + ,^1—G%{E + oj)G^„{E) J dr'Gw(r - r')C^(r' - r)7'^-(q, co; r'). (B4) 

Thus, we find that 

7'^^(q,c^;r)=7'^'^(q,c.)e^'i'- (B5) 

with 

j'^''iq,u) = l+^^^^^^r''{q,^)G%{E + L,)G''^iE)Jdr'CN{r-r')CN{r'~v)e-^^^'^-'''l (B6) 
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FIG. 12: Diagrams for the (scalar) vertex corrections in real space. 




Here we used the identity \/2'itvqTq — (27r^^)A^/4. The integral is equal to 

j dr'CNir ~ r')CN{r' - v)e-^^^^-^'^ = 2^^'"'' [^^ ('z'^'/S)]' 

^ 2^'^oi'lRc), (B7) 

where the second equality holds in the limit of high Landau levels. Neglecting the frequency-dependence and using 
the identities 



4 

2 
1 



GjyGjy — -^ (B8) 






we can solve for 7'^'', and obtain Eqs. (|18|) and p9|) for the vertex corrections. Finally, for finite u we get 

1 

which is used in Eq. (|28|) . 



7"'^(q,^)= . ..2/,w2. „^^M.^ , ..^^w^^ . (BIO) 



APPENDIX C: CORRECTIONS OF ORDER A/uj^ 

In this appendix, we consider the contributions to the triangle vertex to order A/ujc in more detail. To this order, 
vertex corrections of the scalar vertices can be neglected. 

We first consider the case (i) in which both Green functions adjacent to the current vertex are evaluated in Landau 
levels other than N. As mentioned in Sec. lIII(T2l the Green function connecting the scalar vertices should be evaluated 
in the A^th LL up to corrections of order (A/cJc)^. Using the semiclassical expression (|23|l for the matrix elements, 
we then find for the corresponding correction to T^''' the expression 



5rW(q,c.) = ---i-/^(G+-G^) ^ r') \ — l^^^J^_„+l(gi?,)J^_„(gi^, 

ITT 2ttP £mV2 n^N~N+i ^^ ^ (^ ^ "'^'^ {N -n + l)wc 



Lu VN (G^-GJ,)^ 1 /2^ 



27r2^2 ^^^ ^2 Z^k{k + l)\0 



Ettf-tt [M<iRc)Jk+i{qRc) + J-kiqRc)J-ik+i)iqRc)] ■ (CI) 



fc=l 

Using that J-k{z) — (— l)'^Jfe(z), we find that the expression in square brackets vanishes so that (JF^^-' (q, lj) = 0. 
The corresponding contribution to F^"^ (q, uj) takes the form 

<5F('^)(q,^) = ^-l^lniA^' J^„(<zi?,)G+G+ 

1 fi ^^"^ 

= —-^:^^^-gz22'^ki<lRc)iGJ^G^_^ + GJ^GJ^^^). (C2) 

The prime on the sum indicates that only those terms should be kept, in which one of the two Green functions is 
evaluated in a LL different from N. In leading order, 

r+ — -r+ — ^ 
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and hence also the contribution ljC2|) to the triangle vertex vanishes. 

Next, we turn to the contribution (ii) in which the diagrams in Fig. 0] are evaluated to next-to-leading order in 
A/wc while neglecting the vertex corrections on the scalar vertices. Such contributions can arise in particular from 
the self energy entering Gn±i- We first consider the corresponding contribution to Tx ■ (For the purpose of this 
appendix, we choose q || x.) According to the diagrams in Fig. 21 we have for the contribution (ii) 



pW = ^- ^ 



JoiqRc)JiiqRc)2lin{G]^_, - G^+i)[G+] 



+ 12 



TT 27r^2 ^^^"- 
Here, we have already used that to the order under consideration, 

To our order, we then find 



(C3) 



(C4) 



/2N 



TT I 



(C5) 



Comparing with Eq. (|25|) . we find even to order A/uJc that this contribution is cancelled exactly by Tx ■ Thus, there 
is also no contribution of type (ii) to F^; and Tx vanishes to order A/cOc- 

Finally, we consider the contribution of type (ii) to Fj^. Since this is a transverse contribution, we need to consider 
only r^''). In this case, the diagrams in Fig. 21 translate into the expression 



iu} 1 



n 27r^2 ^^^ 



iMqRc)JiiqRcKG+ - G^)2iIm(G^_i + G]^^,)G 



N 



(C6) 



Noting that the leading order cancels from the combination, Gj^_^ + G^y i i, we can simply evaluate Eq. (jC6p for Fj, 
to leading nonvanishing order. This yields Eq. (|38|l for r'^'^/'^'=^ in the main text. 



APPENDIX D: CONTRIBUTIONS TO DRAG FROM DIFFERENT MOMENTUM REGIONS 



We write down explicitly the momentum integrals determining the function I{uj) in Eq. (|67|) . The first integral, 
corresponding to the diffusive range of momenta qRc <S^ 1, 
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(Dl) 



is dominated by the contribution of the "diffusive rectification" , Eq. (|31|l , while the screening is determined by Eq. H34|l . 
The second integral 



Iu{uj) = /ii-i(cj) + /ii_2(w), 



(D2) 



includes the contribution of F^i/'^-) [denoted by /ii_i(w)] and F^'^/'^-) [denoted by /ii_2(w)], Eqs. ^ and l|SS|l . 
respectively, while the screening in In is determined by A''th LL, Eq. (|51|l . 
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The integration domain in the third integral, /iii(ti;), corresponds to the range where the screening acquires its static 
zero-B form, Eq. H53|) . while the triangle vertex is dominated by r('^/"'=\ Eq. H38|l . 
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Let us analyze the first term in In, Eq. (|D3|I . Consider identical layers. The screening is nontrivial and almost 
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vanishes in the vicinity of zeroes Qn of jQ{qRc). The structure of the integral is 



(D6) 



where A ~ lOcj A. We see that the integral is dominated by the momenta close to Qn-, each peak contributing 
~ i?-2gy^A-'^/2, so that the total result 
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is determined by the upper limit where AJ^qRc) ~ 1. 
Similarly, we estimate the second term in In, Eq. (|D4|I . 
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(D8) 



yielding the result of the same order as for Eq. (jD3p . since both integrals are dominated by the upper limit. We note 
that for this reason the same estimate can be obtained by replacing Jq (qRc) , Ji (qRc) by {irqRc)"^. The two terms 
/ii-i and /ii_2 give contributions of the opposite signs to the drag resistivity, since 0{l/qRc) <^ ^\\{B) — r||(— -B), 
while 0{A/ujc) ^ T±{B) = -r±{-B). 
Estimating other terms, we obtain 
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where in the diffusive term R the momentum integration is restricted from below by Qmin = Rc{ijJ / KQaD)^'"^ ~ 
Rc{T / K,QaDY'^ . This infrared cut-off is necessary, since the momentum integral diverges logarithmically at small q in 
the diffusive regime, when Eq. H48|l and Eq. (|34|l are used for the interlayer interaction. The divergence is naturally 
cured when the general formula H47() is employed together with Eq. (|34|l . 

Thus we conclude that at low temperatures T ^ A the total integral is dominated by the contribution of high 
momenta, /m. 
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FIG. 13: Functions Q{x) - Eq. ^, and n{x) - Eq. ((^ZJ, and th e product V'^{x)y\!(x) - Eqs. (jHSJ and (JD^OJ, determining the 
frequency dependence of the "inelastic kernel" /(tj) - Eq. <D15ll . 

resulting in Eq. {TUJ. 

In the case of higher temperatures, A <C T ^ cjc, the main difference is related to the fact that the contribution 
of a single LL to the polarization operator is thermally smeared, yielding an extra factor ^ A/T as compared to the 
second term of Eq. H51|) . as follows from Eq. (|54|l . This changes the upper (lower) limit of integration in /n (/m) 
where A should be replaced by T. Furthermore, in I\\ one should replace Wc/A by lOcJT in the factor related to the 
screening, which is equivalent to multiplying A by A/T in Eqs. IJD6|I and HD8|I . This yields 
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1 
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KT 



(D13) 
(D14) 



We see that for A < T < Wc the contribution of the F^^/^^'^^-term to the momentum integral increases faster than that 
of r^ ''^'='-terin. To evaluate this contribution more accurately, we consider the corresponding momentum integral in 
the whole range of q and include the imaginary part of n(q, cj) into the screening (for simplicity we consider identical 
layers). 
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where V{x) is defined in Eq. (|S5|l and 
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(D15) 
(D16) 

(D17) 
(D18) 



according to Eqs. H54() and (|56|l . The functions Q{x), and 7i(a;) are presented in Fig. 1131 

From the above estimates we know that the momentum integral is determined hy q ^^ lUc/TRc ^ ^/Rc- This holds 
provided A,B ^ uJc/T > 1, i.e. for cosh([£;Ar - ^i\/2T) < (cjc/T)^. On the other hand, uJc/TRc < l/a in the 
ballistic regime. In this case, we can set g^/sinh qa ~ l/a? in Eq. IJD16|) and set the lower integration limit to g = 0. 
Separating the fast and slow variables in Eq. (|D16p . we get [Ji{zn) = 0, Zn — imA- 7r/4] 
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26 
' ^L]\osh'(^^-^]w{u/2A), (D19) 



yv{x) = 



32a2i?2 y^^J \ 2T 

1 r, 2(2^) 



2:27^2(2.) 1^ |x|7^(x) 



n ^ Qix) 

— — arctan ■ 



2 \x\n{x) 



(D20) 



Substituting this result into l|D15|l and integrating the obtained /(w) over frequency according to (|67fl . we arrive at 
Eq. H83|) of the main text. 
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